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The d* -^ drnr partial decay width has been calculated in a wave-function model for d* and d. It is 
found to be smaller than a previous estimate by a factor of 7. A previously proposed dependence 
on dibaryon sizes is confirmed. The large reduction found is caused partly by a change in the d* 
size used and partly by the need to match momenta in pion emissions, a feature not included in the 
previous estimate. 



PACS number(s): 14.20.Pt, 25.40. Ep 

Years ago, Goldman et al. [|l| (referred to below as 
GMSSW) have estimated the dinT partial decay width of 
the dibaryon d* ( J'^ = 3+ , T = 0) to be about 20 eV at the 
dibaryon mass m* 100 MeV above the dmr threshold of 
about 2150 MeV/c^. Though small, this width is of some 
interest as the total pion production cross section of d* 
from deuteron can be related to it. The result estimated 
by GMSSW is 0.1 /xb at an incident pion momentum of 
580 MeV/c. 

To estimate the decay width, GMSSW use an effec- 
tive dd*TnT vertex of the form {g/ml)d*r ^^yd^d'^Trid''TT2, 
where {fj,vX} are symmetrized indices. The effective cou- 
pling constant g can be related to the effective coupling 
constant g-^NA for A -^ Ntt decay by the expression 



g « gl^^{mjMN)^fid* ^ AA)/(d ^ NN) 



(1) 



where the vertex overlap factors are taken to be f{d^ 
NN) - 1 and f{d* -^ A A) - [r* /ra)^ - 2-^, where 
rd (r*) is the deuteron (d*) radius. In this approximate 
treatment, one of the vertex factors can always be set 
to one, since it is only the volume ratio that matters. 
The reduction shown for the second vertex factor repre- 
sents the estimated difficulty of finding the same three 
quarks in each baryon in the smaller size in d* when the 
baryon was originally in the larger bound state that was 
d. This GMSSW approach to the dinT partial width is 
elegant, but it has an Achilles heel, namely the inability 
to improve on the given choice of the vertex factor. It 
is therefore of some interest to determine if the reported 
estimate is reliable. 

Recently, I have estimated the NN and ttNN par- 
tial widths of d* using a perturbative method based on 
a baryon wave function for d* P,p|. The purpose of 
this brief report is to illustrate how this more system- 
atic approach can be applied to the calculation of the 
dirn partial width. The calculated width shows the size 
effect proposed by GMSSW, but depends additionally 
and strongly on the momentum matching required in the 
pion emissions, a feature that is completely missing in 
the GMSSW approach. 

For an S-wave deuteron state described by a realistic 
Bonn C wave function B and a three-term Gaussian ap- 



proximation to a two-center baryon wave function for d* , 
the width calculated here is about 3 eV at m* = 2250 
MeV/c^. This estimate is smaller than that reported by 
GMSSW by a factor of 7. About half of this reduction is 
caused by our use of a smaller volume ratio, while the re- 
maining half of the discrepancy comes from the effect of 
momentum matching among the wave functions required 
in the pion emissions. 

The d* — > diTTT decay in the lowest-order perturbation 
theory can be visualized as the decay of the two off-shell 
A's shown in Fig. n in time-ordered diagrams. The re- 
sulting decay width is, according to the Fermi golden rule 
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where Vi is a pion-emission vertex from a baryon, and 
El and E2 are the pion energies in the d* rest frame. 
G12 is the sum of Green functions for the two diagrams 
labeled a and b. To avoid possible zeros in the energy 
denominators Di, i — a,b, where for example 

Da^m*-Eiiki)-EN{-p-ki)-EAip) , (3) 
I give m* a total width Tfof. Consequently 
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The Green-function factor G^2 depends on the an- 
gles. However, if it is approximated by a suitable angle- 
averaged value (0^2), and if the baryons in both initial 
and final dibaryons are in relative orbital S states, the 
angle integrations shown in Eq. ^ can be performed 
analytically to give 
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Although the baryon-baryon relative wave functions 
of both d* and d will eventually be expressed as sums 



of Gaussians, it is sufficient to show the result for the 
spin- and angle-averaged squared matrix element in the 
integrand for single Gaussians such as 
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for d* . The result obtained after the integration over the 
internal momentum p shown in Fig. n^ is 
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where 
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Here Ta = 120 MeV is the A width, k*{E^) is the mo- 
mentum of the decay pion (recoiling nuclcon) in the rest 
frame of the decaying A, and 



ao = r^/3 = 0.12 fm^ 



(9) 



is the parameter, taken to be the same in N and A, de- 
scribing the Gaussian wave functions of quarks inside 
these baryons. The baryon form factors in the pion- 
emission vertex give rise to the ao term in a. 

When the baryon wave function of the dibaryon d* 
is approximated by a single Gaussian, the parameters 
f3* and the m.s. baryon momentum < p^ > inside the 
dibaryon are related to the dibaryon radius r* as 



p* = y^/r* , < p2 > = 
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The average baryon energies to be used in the Green- 
function factor < GI2 > are 



i^A(p) 



Ml+ < p2 > 



Af2,+ < p2 > +A:2 . 
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Certain features in our results are worth pointing out: 
The /3-dependent factor in Eq. (0) can be expressed ap- 
proximately in terms of baryon radii as 



S2 



Td 
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This is the only feature included in the GMSSW ver- 
tex factor. Note however that if we use r^ = 2.0 fm 
and r* = 0.7 fm, then this volume ratio is smaller than 
that used by GMSSW by a factor of 3. It is possible 
that GMSSW have used a larger volume ratio in the 
expectation that in the quark-delocalization model the 



quarks from a baryon constituent of the dibaryon are 
more spread out over the dibaryon than in a conven- 
tional baryon-baryon bound state. Unfortunately, this 
expectation has not been quantified in GMSSW. 

A quantitative study of this effective size effect of de- 
localization on the dTTTT decay width is not easy because 
delocalizaed quark wave functions IIM] are so compli- 
cated when antisymmetrization and angular momentum 
projections are made. In addition, delocalized short- 
distance wave functions will also have to be used in the 
deuteron for overall consistency. In the more literal ap- 
proach taken here, this important question will be left as 
an open problem for future study. On the other hand, 
certain other features of the problem, such as an increase 
in the dibaryon size, can be studied rather easily in the 
simple model used in this paper. 

The remaining features of Eq. ([7|) are not included in 
the GMSSW vertex factor. They describe the require- 
ments of momentum matchings on pion emissions. Note 
that the parameter B^ is dominated by /3*2, which is pro- 
portional to the m.s. baryon momentum contained in d* . 
For r* = 0.7 fm, B^ = 0.29 fm2, which is larger than ao 
by a factor 2.4. This shows that the dependence of the 
parameter a on baryon sizes in the inelastic baryon form 
factors is significantly weaker than the dependence on r* . 
If the A radius ta is different from the nucleon radius r^, 
the parameter ao should be replaced by 
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Since ta is expected to be larger than r^ by only 10 
(~ 5) % in the Massachusetts Institute of Technology 
bag model (many potential models |^), the overall 
effect in the parameter a is only 3(2)%, which is quite 
negligible. In contrast, results will be shown below where 
the parameter a changes by a factor of almost 2. I there- 
fore conclude that the use of the same baryon size for 
both A and N is justified both here and in the "delo- 
calization" model of Mpl. Without this simplification 
of equal baryon sizes, the delocalization model would be 
even harder to execute technically. 

The total width Ttot of d* is also needed in the calcu- 
lation. It is known to increase with increasing d* mass 
m*, rising from about 1 MeV at to* = 2100 MeV/c2 to 
about 10 MeV at to* = 2350 MeV/c2 |J|. However, 
the presence of Ttot is unimportant at the low end of 
this mass range, because the energy denominators Di are 
large there. Beyond the AA^tt threshold at 2310 MeV, 
however, the energy denominators could become small; 
the choice of Ttot then becomes important. For this rea- 
son, I use a nonzero width in the calculation, but choose 
for simplicity a single constant value of Ttot = 10 MeV 
appropriate to the high end of the mass range studied 
here. The effect of using smaller values of Ttot will be 
explicitly shown below. 



The single S-state Gaussian wave function for the 
deuteron depends on the parameter /3 = y/3/8/rd, where 
rd = 1.967 fm is the radius of the deuteron wave func- 
tion. The calculated dTTTr partial width for r* = 0.7 fm 
is shown in Fig. as a thin solid curve for a range of 
m*. The value at m* = 2250 MeV/c^ is 2.0 eV, ten 
times smaller than the estimate given by GMSSW. Be- 
sides the decrease by a factor of 3 already pointed out 
and discussed previously, there is a remaining reduction 
by another factor of 3 which should be attributed to the 
technical improvements made in the present calculation. 

For a more realistic S-state wave function, I use the 
three-term Gaussian fit to the Bonn C deuteron S-wave 
wave function obtained in g. The fitted S-state prob- 
ability of 94.34% (versus 94.39% for the original Bonn 
C wave function) has been renormalized back to 100% 
for a pure S state. The resulting dnn partial width is 
shown in Fig. || as a solid (dashed, long-dashed) curve 
forr* =0.7(0.9,0.5) fm. 

These curves show that the decay width for the real- 
istic Bonn C wave function is larger than that for the 
single Gaussian wave function rather uniformly over the 
mass range. At m* = 2250 MeV, the increase is by a 
factor of 1.9 at r* = 0.7 fm. Since the deuteron radius 
is essentially the same in these calculations, the change 
must have come entirely from the increase in the high- 
momentum components in the Bonn C wave function. 

Figure also shows results calculated with the same 
Bonn C wave function but different d* radii. In the sim- 
ple volume scaling model of GMSSW, the decay width 
for r* = 0.7 fm would be increased (decreased) a factor 
of 2.2 (2.7) when r^ is increased to 0.9 fm (decreased to 
0.5 fm). The actual calculated factor at m* = 2250 MeV 
turns out to be 1.6 (2.7) for the Bonn C deuteron, and 1.9 
(3.3) for the single-Gaussian deuteron. Thus the overall 
volume scaling effect proposed by GMSSW appears to 
be present, but in a more complicated form that also de- 
pends on other details of the wave function. The main 
reason for the complication is that the parameter a also 
changes drastically with r* — from 0.40 fm^ for r* = 0.7 
fm to 0.28 (0.58) fm^ for r* = 0.5 (0.9) fm. 

The additional sensitivity to the high-momentum com- 
ponents of the deuteron wave function shown in the 
present calculation suggests that the deuteron D-state 
component, though constituting only 5% of the deuteron, 
might have a disproportionate effect on these decay 
widths. Its inclusion is easy to visualize, but tedious 
to execute. I believe that its inclusion will not change 
the order of magnitude of the partial widths estimated 
here — after all, the effect of changing the S-state wave 
function from the smooth single Gaussian to the realistic 
Bonn C wave function is an increase by "only" a factor 
of two. 

I now turn to the effect of using smaller total decay 
widths in the calculation. The results for Ttot = 5 (2.5) 
MeV are shown in Fig. as a dashed (dotted) curve. 



and compared with the corresponding result for Ttot — 
10 MeV from Fig. ||, reproduced here as a solid curve. 
The calculated width can be seen to develop significant 
dependence on the total decay width only well above the 
ANn threshold of 2310 MeV. 

The single Gaussian wave function for d* is of course 
a rather crude approximation to the model of "delocal- 
ized" quarks in d* wH . The least improvement one could 
make is to use a two-center Gaussian wave function for 
each baryon proportional to exp[-(x - s)^/2] -I- exp[-(x 
-I- s)^/2], where x is the (dimensionless) relative baryon- 
baryon coordinate and 2s is the separation between the 
two centers. The projected S-state component of this 
relative wave function has the form 



ipoix,s) = Nq 



1 r 



-(x^sf/2 



-{x+sf/2 



(14) 



The lowest-order effect of s can be included in the single 
Gaussian wave function by matching the m.s. radius 
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For small s, one sees the separate contributions from the 
zero-point motion of the baryons in the two harmonic 
oscillator potentials and from the separation 2s of these 
potentials. If one starts with quark wave functions de- 
scribed by a size constant b = 0.6 fm, the relative baryon- 
baryon motion for the six-quark state will be described 
by the size constant br = b/V6. Hence r* — br (x^)^'^. 
For a potential separation 25* = 1.40 fm calculated for 
the d* @,§, s = S/br = 2.86, and hence r* = 0.72 fm. 
This is close to the middle of the three values used in Fig. 

i . 

It is easy to go beyond this rough approximation and 

actually fit the momentum wave function of the projected 
two-center Gaussian, namely ^ exp[— (6rp)^/2]j'o(pS'), as 
a sum of three Gaussians 

3 

'^p2c{p)^'^c^■^pi{p), (16) 

where ipi{p) is a normalized Gaussian wave function. In 
order to emphasize the stronger high-momentum com- 
ponents introduced by the oscillating factor joipS), the 
range parameters are obtained by minimizing the per- 
centage m.s. deviation. The fitted result is 



7 -(71,72,73) =0.9961(1, 1.25, 1.46)/62 



(17) 



where 7i = 2/3^^. The expansion coefficients, renormalized 
from the fitted normalization of 1.0019 back to 1, are 

c = (ci,C2,C3) = (14.8311,-27.4124,13.3505) . (18) 

The drnr partial decay width can now be calculated with 
this improved d* wave function but with everything else 



treated in the same way as the single-Gaussian case with 
r* = 0.7 fm. The results, shown as solid circles in Fig. 
g, are about 20% smaller than those for the single Gaus- 
sian d* wave function. The radius of the fitted d* wave 
function, at 0.72 fm, is actually marginally larger, but 
the decrease in the calculated decay width shows that it 
is the additional high-momentum components in the d* 
wave function that dominate the result. 

Although the present calculation is a significant im- 
provement over that of GMSSW, it is also not a quan- 
titative calculation. Too many approximations have 
to be made to reduce the problem into manageable 
form. The most important features that should be in- 
cluded in a realistic calculation include the following: (1) 
higher-order Feynman diagrams, especially those involv- 
ing intermediate- and final-state interactions; (2) quark 
antisymmetrization between baryons and quark delocal- 
ization effects in both d* and d; (3) deuteron D state and 
perhaps better short-distance wave functions for both d 
and d* , in both NN and exotic channels; (4) improved 
treatment of the energy denominator, requiring the use of 
full angle integrations; (5) better treatment of the AA^tt 
vertex; (6) a better choice of the full decay width, but 
only when the dibaryon mass is above the ATVtt thresh- 
old. It is clear, however, that these improvements will 
require much more extensive calculations than those un- 
dertaken here. 



Figure Captions 

FIG. 1. The two leading-order time-ordered diagrams for 
the decay d* — » dim. 



FIG. 2. The d* —* dim partial decay width rd,r7r as a 
function of the d* mass m* for single-Gaussian and Bonn 
C deuteron S-state wave functions, and for single-Gaussian 
and two-center d* wave functions. The number shown in the 
legend is the r* radius in fm. 



FIG. 3. The d* — > div-R partial decay width Vd-mr as a 
function of the d* mass m* for Bonn C deuteron S-state wave 
function and single-Gaussian d* wave function (with r* — 0.7 
fm) when the input total width of d* is varied. 
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